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1 Introduction 


The Einstein-Maxwell theory forms the basis for other gravitational-electromagnetic 
theories. The Einstein-Maxwell theory arises from the Einstein-Hilbert gravitational 
action plus the Maxwell action. It is minimally coupled because there is no coupling 
in the Lagrangian between the Maxwell part and the curvature part. It also gives 
equations which are second order in the derivatives of the metric (as opposed to 
higher order), because the Lagrangian does not contain generic products of curvature 
terms (the second derivatives of the metric that might appear in the Einstein-Hilbert 
gravitational Lagrangian form a divergence of some vector and do not contribute to 
give equations of higher order). In addition the Einstein-Maxwell theory is linear 
in the electrodynamics, which means that the Maxwell Lagrangian is quadratic 
in the Maxwell tensor. The property that most interest us here is the coupling 
between the electromagnetic and the gravitational parts. Thus we are led to classify 
gravitational-electromagnetic theories in a useful way into two classes, according to 
the theory is minimally coupled or non-minimally coupled. 

The first class is minimally coupled gravitational-electromagnetism. It can have 
different subclasses. One can subdivide into two subclasses, whether the correspond¬ 
ing electrodynamics is linear or non-linear. One can then also subdivide into new 
sublacsses according to the gravitational action, whether the gravitational part yields 
a second order theory (such as Einstein-Hilbert theory), or a higher order theory. For 
instance, a minimally coupled theory, linear in electrodynamics, and second order in 
the gravity part is the standard Einstein-Maxwell theory [T]. There are many exact 
results in the framework of this theory, such as the Reissner-Nordstrom solution for 
a charged black hole, gravity wave solutions in electrovacuum, cosmological models 
with a magnetic field, to name a few. Another instance is a minimally coupled 
theory, with non-linear electrodynamics, and second order Einstein gravity. The 
well-known models of Born and Infold [ 2 ] and of Heisenberg and Euler 0 when cou¬ 
pled to gravity belong to this subclass. One of the most interesting problems in this 
theory is the search for the regular, non-singular, black holes. This search, started 
by Bardeen in 1968 0, and developed by many authors (see, e.g., 001311]) led 
to the recent success of Ayon-Beato and Garcia umm in finding exact solutions 
to the Einstein equations coupled with specific four-parameter models of non-linear 
electrodynamics. And of course there are the other cases of higher order theories 
coupled to Maxwell, or to non-linear electrodynamic theories. 

The second class is non-minimal coupled gravitational-electromagnetism. It can 
be subdivided according to whether the corresponding electrodynamics is linear 
or non-linear. Now, in non-minimal coupled models one can no longer divide, a 
priori, into second order and higher order theories, since by definition curvature 
terms appear in these models which in principle give rise in general to higher or¬ 
der terms in the equations. This second class includes non-minimal equations for 
the electrodynamics containing couplings with the Riemann and Ricci tensors and 
the Ricci scalar. This class is very wide and comprises several subclasses, such 
as: non-minimal linear electrodynamics plus Einstein-Hilbert term, non-minimal 
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non-linear electrodynamics plus Einstein-Hilbert term, non-minimal linear electro¬ 
dynamics plus Einstein-Hilbert and other pure curvature terms, non-minimal non¬ 
linear electrodynamics plus Einstein-Hilbert and other pure curvature terms, and 
others. All these subclasses of models belong to this second class since they have 
one specific feature: the Lagrangian contains an interaction part with specific cross¬ 
terms, including scalar products of the Riemann tensor and its convolutions, with 
the Maxwell tensor. 

Our goal is to study this last second class. This class is of great interest, since the 
appearance of cross-terms in the Lagrangian leads to modifications of the coefficients 
involving the second-order derivatives both in the Maxwell and Einstein equations. 
This means, in particular, that gravitational waves can propagate with a velocity 
different from the velocity of light in vacuum, in a similar fashion as electromagnetic 
waves propagate in a material medium. This new added feature has many inter¬ 
esting applications in various systems and models, such as cosmological scenarios, 
gravitational waves interacting with electromagnetic fields, and charged black holes. 
More specifically, in cosmology the evolution of the gravitational perturbations may 
have another rate and scale. In astrophysics, the interaction of gravitational with 
electromagnetic waves may lead to time delays in the arrival of those waves, and the 
gravitational waves themselves would change their own properties in a form notice¬ 
able in gravitational wave detection. It also leads to important modifications of the 
electromagnetic and gravitational structure of a charged black hole. 

First we consider the simpler case of minimal coupling in the electromagnetic 
and gravitational parts, in Section 2. Then in Section 3 we study non-minimal 
coupled theories. In section 3.1 we obtain the structure of the master equations of 
the non-minimal gravitational-electromagnetic theory, for both non-linear and linear 
electrodynamics. In section 3.2 we consider in detail the linear version of the theory. 
In section 3.3 we briefly study an example. 

2 Minimal coupling of gravity and electromag¬ 
netism 

2.1 General formalism 

In order to explain the novelty of our approach, let us first introduce the nomen¬ 
clature in the well known case of gravitational-electromagnetic theories minimally 
coupled. We will consider, generically, high-order terms in the gravitational part, 
and non-linear terms in the electromagnetic part. The action functional is jjfj 


S = j d A xy/^g C 


min ? 


(i) 


where, 
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g is the determinant of the metric tensor g ik , and £ m i n is the Lagrangian for the 
minimally coupled theory. It is composed of two distinct parts, which do not cross, 
the Lagrangian £, related to the metric field, and the Lagrangian £, related to 
the electromagnetic field, 
geometrical scalars only, 


The Lagrangian £ R ik R lk , R ikmn R 


jikmn 


contains 


p p pifc p pifcmn /o'. 

i ikmn-Li' 7 ••• j ('-'J 

where R is the Ricci scalar, Ri k is the Ricci tensor, and Rikmn is the Riemann tensor. 
The constant k is equal to k = 

The Lagrangian £(/(n), ip 2 )) is an arbitrary function of the quantities ipi) and 
/p 2 ). /(n) and /( 12 ) form a first set (first subscript 1) of electromagnetic held invari¬ 
ants. This first set is composed of two invariants (denoted in the second subscript), 
the first ipq and the second ip 2 ) invariants. These invariants are quadratic in the 
anti-symmetric Maxwell tensor F lkl and given by 

i m F-F it F‘\ (4) 

The asterisk denotes the dualization procedure, defined as follows 

F *ik = }_ e Ms Fis _ ( 5 ) 


Here e lkls = —£= e lkls is the Levi-Civita tensor and e tkls is the completely anti¬ 
symmetric symbol with e 0123 = — eoi 23 = 1- The Maxwell tensor satisfies the condi¬ 
tion 

V k F* ik = 0, (6) 

where V& is the covariant derivative. Equation can also be written as WiF k i T 
ViF ik + V k F h = 0. Due to ®, the Maxwell tensor may be represented in terms of 
a four-vector potential Ai as 


Rk = — V k Ai 


dA k _ dAj 
dx 1 dx k 


(7) 


Now, the variation of the action functional m with respect to the four-vector po¬ 
tential Ai gives the minimal vacuum Maxwell equations 


Vfc 


F ik , d£ 
dl{ 11 ) dl { 12 ) 


jj*ik 


0 . 


( 8 ) 


On the other hand, the variation of the action functional m with respect to g lk 
yields the gravitational equations 


Ein^ T ik , 

K 


(9) 
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where Ein^ is the corresponding non-linear generalization of the Einstein tensor, 
G ik = Rik ~ \Rgik ■ The tensor T ik , defined by 

1 8 C 1 r)C 

Tik = A - ^—F in F k n - ^^—(F*F k l + F u F * k l ), ( 10 ) 

Z G'-qn) ZOl( 12) 

is the symmetric stress-energy tensor of the electromagnetic held in vacuum. The 
tensor T\ k is conserved in accordance with the Bianchi identities 

V k T ik = 0 . ( 11 ) 

2.2 Example: linear Einstein-Maxwell theory 

Before we leave this section, we give in this subsection, as a simple example, the 
usual linear Einstein-Maxwell theory. It can be obtained from equations m-m 
when the gravitational Lagrangian is given by the Einstein-Hilbert term — only, and 
£(J(iu, Ifi2)) = T(n) . The relations (IHll- lTUl) reduce, respectively, to 



O 

'&H 

> 

(12) 


^ ri _ t^O) 

^ ik > 

K 

( 13 ) 

and 




rp(0) — 1 77> Tpmn T? r?-n 

ik — ^9ik-r mn^ inF k 

( 14 ) 

where the superscript 

(0) denotes that the tensor 

is the simplest part of a 


more general electromagnetic stress-energy tensor. Such a formalism describes a 
minimal coupling of gravitation and electromagnetism, since the right-hand-side of 
the Einstein equations ©, as well as the Maxwell equations © and © contain 
metric couplings and covariant derivatives only, while the curvature tensor appears 
exclusively in the left-hand-side of ©• 

In the following section, and in contrast to the minimal gravitational-electromag¬ 
netic equations discussed in this section, we consider in some detail, along the sec¬ 
tion, a non-minimal gravitational-electromagnetic theory both with non-linear and 
linear electrodynamics, generalizing the Einstein-Maxwell theory and other minimal 
theories. This approach deals with self-consistent modifications to both the Einstein 
and the Maxwell equations. 
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3 Non-minimal extensions 

of the Einstein-Maxwell Lagrangian 

3.1 Full formalism and equations 

3.1.1 Invariants containing the Maxwell and dual-Maxwell tensors cou¬ 
pled with the Ricci scalar and the Riemann and Ricci tensors 

Let us introduce the invariant scalars, quadratic in the tensors F lk and F* k and 
containing the Riemann and Ricci tensors and the Ricci scalar. These scalars yield 
cross-terms, and are the appropriate quantities for the description of non-minimal 
interactions. They can be formally divided into five sets. The first set is the trivial 
one, containing 7(n) and /( 12 ) alone, as described before. The second set contains 
/(ip and /( 12 ) multiplied by R , 

ini) = jg ,m g k "F< k F mn , 7,22, = jg' m g ln F; k F mn . (15) 

The third set includes the Ricci tensor R mn , 

I(3i) = l -R vm g kn F lk F mn , J (32) = l -R im g kn F: k F mn . (16) 

The fourth set is based on the convolutions of the quadratic combinations of F lk 
and F* k with the Riemann tensor 


R 41 ) 


= -R 


• ikmnp p 


ik”mn ? 


t — L TDikmn 77 * 77 

7 ( 42 ) = 7; K ^ik^mn 


(17) 


The invariants I( 21 )—1(42) are chosen to be linear in the curvature. Note also that the 
scalar R rrn g kn F* k F* nn can be reduced to a linear combination of 7( 2 p and 7( 31 ), and 
the scalar \R lkmn F* k F^ n can be represented as a linear combination of I( 21 ), 7( 31 ), 
7(2H). Finally, the fifth set includes the various scalars nonlinear in the curvature. 
Below we introduce a few of them, 


I (51) = -g im g kn F lk F mn f R , 

1 

T — Dim. pfcn 771 771 

4 ( 53 ) — 2 44 F i k r mn i 

t _ ^ ryikab p l. ' rrnmn 

4 ( 55 ) — 2 ^ -F^abrnn-F ik-F > 

t _ ^ rg*ikab p* I' jgmn 

4 ( 57 ) = -p K abmn r i k r , 

t _ ^ r->ikab p 

4 ( 59 ) — 2 44 Fiabcd 


1 

T — rd m r, kn T?* T? T? 

4(52) — 2 ^ 9 FfaF mn r R , 

1 

T — nim pfcn p* p 

4(54) — 2^4 44 F ik r mn , 

7 _ ^ rjikab p 771* 77177177 

4 ( 56 ) — 2^4 F^abmnF i k F , 

7 _ 1 T}*ikab p* p* P 7717 ! 

4(58) 2^4 Fi'abmn, ik, > 

Tjcdmn p p 
44 F{ k F mn , ... . 


(18) 

(19) 

( 20 ) 
( 21 ) 
( 22 ) 
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In equation (fTHlh /r and Fr denote arbitrary functions of the all possible indepen¬ 
dent non-linear scalar invariants of the gravitational held, such as R 2 * 4 , R mn R mn : 

td jDikmn td* TDikmn 
■R'ikmn**' i ri ikmn ri > ••• 

Thus, the non-minimal Lagrangian can be written in the form 


r . =£ 

^nnn-rmn 


R 


K 


"? Rn 



(23) 

This non-minimal Lagrangian is 17(1) gauge invariant since it contains the Maxwell 
tensor F lk only, and does not include the potential four-vector A\ 


3.1.2 Non-minimal non-linear electrodynamics 


The variation, with respect to the 4-vector A k , of the action functional with La¬ 
grangian (El yields the equation for the non-minimal non-linear electromagnetic 
held, 

= 0 , 

where H lk is the induction tensor given by 


Tjik \ likmn r? \ -\\nkmn r£* 

ti v r rnn ~r Vv -t rrii 


(24) 

(25) 


with 


V 


ikmn _ 


_ 2_ ^gimgkn gkmgin\ 


dc d£ „ 

+ -—R 


dL 


( 11 ) 


dL 


( 21 ) 


+ 


1 B£ 

( rfim „kn ryin km , nkn im rykm in\ 

-7 \R 9 -R 9 +R 9 -R 9 J 77 — 

4 01 ( 31 ) 


+R 


ikmn 


9£ 


dL 


(41) 


+ ... , 


and 


■yyikmn _ 

1 


_ ^gimgkn gkm gin^ 


9£ 9£ „ 

+ ttt —R 


dL 


( 12 ) 91(22) 

im Azn r>in „km 


+ 


(26) 


t~)( „im „kn „km An\ ( TDim „kn r)in „km . r>kn Am r>km An\ 

ti{g 9 -9 9 ) ~ \ti 9 -tig + ti g — K g 


9£ 


J dL 


(32) 


+ 

+ ... . 


73 1 

r>ikmn , ( im kn „km in\ 1 ( 737771 „kn r>in „km . r>kn Am rykm. in\ 

R + ~£\9 9~9 9 ) - -^yR 9 ~ K 9 +K 9 -R 9 ) 


9£ 

dl(i2) 

(27) 


3.1.3 The generalized higher order Einstein equations 

The variation, with respect to the metric coefficients g lk , of the action functional 
with Lagrangian El yields the mon-minimal extension of the Einstein equations, 


1 


0 = Em ik --£g ik + 

2 \dL( ip 


9 £ „ 9 £ 

+ R- 


9L 


FinF-J 1 + 


( 21 ), 
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1 ( 5 £ 


2 \()L 


+ R 


5 £ 


( 12 ) 


5 / 


( 22 )/ 


( BC BC 

(FuF k l + F^F?) + R ik /(ii)^— + /(i2)T^— 

\ Ol (21) Bl(22) , 


l dC, dC 

— V,V fc ) | /(li)^y-b /( 12 ) 


( 21 ) 


< 9 /, 


+ 


( 22 ), 


1 r F ,n (i2 i ,F fcn +i2 w F in )+i2 mn F im F fcn ] +^ ifc V m V JJ^F mn Fl 


2 5/(31) 

-|V m V m 

^v m |v' 

>‘{ v - 

1 5 £ 

(32) 


5 £ 

5/(31) 
5£ 


FinFjJ 1 




.5/(32) 

5 £ 


(f;;f-+z^ 


=m^ 


5 / 


(32) 


(F fc * n F to +F fcn F* ta ) 


/ _ 5 £_ 

V 5/(31) 
+5 , ifcV/ 

+Vfc 


I +V, 


, 5/(31) 

/ 5 £ 


5 / 


■F F 

7.77. -*■ 


In 


(31) 


5 £ 


5/(32) 

5 £ 


^p*mn pi | pmnp*l ^ 


5 / ( 


(32) 


(F* n F ln +F in F 


*ln 


+ 


16 5 / 


[F(F;;F- + F kn F: n + F in F k n + F fen FT)+ 

"b ( F£ n Rkm + F kn R im ) F mn + ( F in R km + F kn R im )F 

' 5 £ 


Q O /• -1 

f /s (ft£W + + -V m v 

4 02 ( 41 ) 2 

3 5 £ 


5 / ( 


(41) 


(F; n Fr + FT^T) 


+ 


( 7T'*m p ZT’fcS J Z7'*m p pts I T7'' m p p*/s . p-m p t?*Ls\ . 

ttkmls-T + -T/c riimls-r ~r Pi Ftkrnls-F + rCimls-F ) + 


8 5/(42) 

1 5£ 1 n 

"^9ik— -1(42) + T V m V r 


5 / 


(42) 


5£ 


5/, 


(42) 


/^*ra^-n I ^*ra^7>-n _i_ p-m p*n i p-m p*n\ 
vz k k % % k k % J 


( 28 ) 


These equations can be rewritten in the well-known form Ein^ = but it 

is not the canonic form when one is dealing with general non-minimal non-linear 
electrodynamics. The reason is the following: even if the Lagrangian for the pure 
gravitational field is of the Einstein-Hilbert form, the equation (OKI) contains higher 
order derivatives of the metric, coming from the curvature tensor in terms containing 
non-minimal scalars, J! c . Thus, a generic non-minimal non-linear electrodynamics 

^(a 6 ) 

is associated with a higher order gravitation. One can then ask the question, whether 
or not non-minimal non-linear electrodynamics models exist for which gravity is of 
second order. We believe that this is possible for a special choice of the dependence 
£(/(aft)) and for specific symmetric space-times. Below we consider a simple example 
confirming such idea. 
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3.2 Non-minimal coupling models, with the coupling linear 
in the curvature, in Einstein-Hilbert gravity 

3.2.1 The action 

A special case worth of discussion is when one restricts the above theory to a La- 
grangian that is Einstein-Hilbert in the gravity term, quadratic in the Maxwell 
tensor and the coupling between the electromagnetism and the metric is linear in 
the curvature. Thus the theory may contain the invariants I(u), Ipi), I (31), /(41), 
only. Such a Lagrangian takes the form 

JD I 1 

r = _ 4 - _ F F mn + - Y ikmri F i F 1291 

**' ' mn- 1 1 0 A ± ik ± mn ? ) 

K Z Z 

where the quantity x lkmn is the susceptibility tensor. The origin of such a terminol¬ 
ogy is the following. One obtains from the Lagrangian m with the definition (ESI) 
that the induction tensor H lk and the Maxwell tensor F mn are linked by the linear 
constitutive law (see, e.g., Baua) 

H ik = pik + x ikmnp mn _ ( 30 ) 


Another important tensor, appearing in the electrodynamics of continuum media, 
is the polarization-magnetization tensor M lk , defined by 


4ir M ik = H ik - F ik , 


(31) 


and equal to 


4tt M ik 


_ ikmn p 

A r mn 1 


(32) 


according to In the standard terminology of continuum electrodynamics na 

m the proportionality coefficients x lkmn form the so-called susceptibility tensor. 
Generally, it has the same symmetry of the indices transposition as the Riemann 
tensor, and has 21 independent components. In our case the susceptibility tensor is 
linear in the curvature. 


3.2.2 Susceptibility tensor 

According to the specifications above, the susceptibility tensor has to be of the form 


X ,tmn = 4^(9 


im „kn An „km\ . ( T)im „kn TDin „km \ rykn Am T)km An\ i -nikmn 

9 ~9 9 j+yl^ 9 9 g —K g j+q^n 


( 33 ) 

The parameters q\, q 2 , and q 3 are in general arbitrary. They have to be chosen 
by some ad hoc constraint, phenomenological or otherwise. For instance, the La¬ 
grangian of the type given by equations (EE)1) and (1H31) . with q t = q 2 = 0, q 3 = —A, 
and A a constant, has been proposed phenomenologically by Prasanna in the con¬ 
text of non-minimal modifications of the electrodynamics ma ng. Some general 
phenomenological properties of the Lagrangian El and El have been discussed 
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by Goenner in (T^j. The problem of a phenomenological introduction of non-minimal 
terms into the electro dynamic equations has been exhaustively studied by Held and 
Obukhov PS]. Drummond and Hathrcll EDI have made a qualitatively new step, 
they obtained modified Maxwell equations from one-loop corrections of quantum 
electrodynamics in curved spacetime. Their model is not phenomenological and 
corresponds to the Lagrangian (GHD and (El with specific choices for q\ , q 2 , and 
q 3 , which involve the fine structure constant and the Compton wavelength of the 
electron. A quantum electrodynamics motivation for the use of generalized Maxwell 
equations can also be found, for instance, in the work of Kostelecky and Mewes 
EH- Accioly, Azeredo, Aragao, and Mukai m used the Prasanna electrodynamic 
equations to construct a special example of a conserved non-minimal effective stress- 
energy tensor containing the Riemann tensor. Exact solutions of master equations 
of non-minimal electrodynamics in a non-linear gravitational wave background were 
obtained and discussed in m-m- 

The susceptibility tensor x lkmn has the same index symmetries as the Riemann 
tensor R lkmn , Its convolutions yield 

g k nX ikmn = R im (q 2 + q 3 ) + \Rg irn (3qi + q 2 ), 
gkngimX lkmn = R(6qi + 3g 2 + q-.i ). (34) 

The coefficients q \, q 2 , and q 3 are considered to be independent phenomenologi¬ 
cal parameters. They introduce specific cross-terms, which describe non-minimal 
interactions of the electromagnetic and gravitational fields. Thus, one has a three- 
parametric family of non-minimal models. We now consider three specihc variants 
in the choice of the set q\ , q 2 and q 3 , and see how it influences the expression for the 
susceptibility tensor. 


(a) The susceptibility tensor is proportional to the double-dual Riemann tensor 

The gravitational analogue of the dual Maxwell tensor F* k , is given by the double¬ 
dual Riemann tensor 


> abed , 


'A ikmn — -**'ikmn — ikabedmn • 


(35) 


The analogy is due to the similarity of the identity \7 n F* n = 0 for the Maxwell 
tensor, with the identity V n Gikmn = 0 for the double-dual Riemann tensor. The 
convolution of the double-dual Riemann tensor gives the Einstein tensor 


l kn Gikmn = Ri 


7,Rgi 


(36) 


Now, the double-dual Riemann tensor is given by 


g 


ikmn _ 


TD 

_ f im kn 

2 yg g 


An „km\ . ( T>im Azn T)in „km . r>kn Am 

9 9 J + 9 9 9 


km An\ 


R g 


-R 


•ikmn 


. (37) 
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Thus, if one imposes that the susceptibility tensor is proportional to the 

double-dnal Riemann tensor, i.e., 

x ikmn = q gikmn , ( 38 ) 

one obtains from equation © a one-parameter model with the following values for 
qi, q 2 j and q\ — q-j, — —q, and q 2 = 2 q. This can also be written as, 

qi + Q 2 + <? 3 = 0 , 2qx + q 2 = 0 . (39) 

For this one-parameter model the non-minimal Lagrangian da) can be rewritten 
in terms of the Ricci scalar, the Maxwell tensor, the dual Maxwell tensor and the 
standard Riemann tensor as follows, 

c = | + \F mn F™ + | R itmn F* t F^ • (40) 


(b) The susceptibility tensor is proportional to the Weyl conformal tensor 
The Weyl tensor is given by 

E> 1 

C ikmn T)ikmn . { Am „kn in „km\ 1 ( r>im „kn r>in „km . r>kn Am r>km An\ ( aa\ 

= ft + ~^{9 9 -g 9 )~ 2 ^ K g ~ K g +K g ~ K g v 4i J 

It has vanishing trace, i.e., gk n C lkmn = 0. If one imposes that the susceptibility 
tensor yf kmn is proportional to the Weyl tensor, i.e., 


ikmn _ nikmn 

A — 9 U J 

one obtains from equation m that 

3qi + q2 — 0 , ?2 d ?3 = 0 . 


(42) 


(43) 


This is also a one-parameter model for which one can easily explicitly give the non- 
minimal Lagrangian (ED- 


(c) The susceptibility tensor is equal to the Drummond-Hathrell tensor 

Drummond and Hathrell 1201 have obtained modified Maxwell equations from 
one-loop corrections in quantum electrodynamics in curved spacetime. Their model 
corresponds to the Lagrangian EHD,® with the following coefficients 

a\ 2 

2<?i-g 3 = 0, I3q 1 + q 2 = 0, q x = - —- , (44) 

1807T 

where a is the fine structure constant and A e is the Compton wavelength of the 
electron. This is also a one-parameter model for which one can easily explicitly give 
the non-minimal Lagrangian (HHl- 
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3.2.3 Non-minimal constitutive equations for the electromagnetic field 

The relation 4HD is of the type of a linear constitutive equation 1121 fUi'j 

H ik = C ikmn F mn , (45) 


where the material tensor C lkmn links the induction tensor with the Maxwell tensor. 
Comparing m with (USD one finds 


(jikmn _ 


im „kn 


;(g im g 


in km\ _i_ ikmn 

y y ) ' X. 


(46) 


The material tensor C lkmn describes the properties of the linear response of the 
material to an electromagnetic field, and contains the information about dielectric 
and magnetic permeabilities, as well as about the magneto-electric coefficients m 

m . 

Using the medium four-velocity U l , normalized such that U l Ui = 1, one can 
decompose C lkmn uniquely as 


C 


ikmn 


1 ^imjjkjjn _ £ injjkjjm + £ knjjijjm _ £ kmjjijjn j + 


+ 7 


-f (/i )i s r s + 7f w (U”V z n - f/ n iy m ) + - U k u t \) . (47) 


Here e im is the dielectric tensor, (p 1 ) pg is the magnetic permeability tensor, and 
u p m is the magneto-electric coefficients tensor. These quantities are defined through 


un 


V, 


pq 

m 


P ■ 


2C ikmn U k U n 
-\g P ikC ikmn ri mn q , 

„ /-iikmnjj _ t t imkln „ 

Vpik ^ Un Uk^ Vlnp • 


(48) 


The tensors g mn i and rf kl are anti-symmetric tensors orthogonal to U l and defined 
as 

Vmni = e mnls U s , rf kl = e ikls U s . (49) 

They are connected by the useful identity 


7 Jkp v _ VkL TJTJS A 4 A k _ A. A k 

'I 1 Imnp u mns u l u L -^n L -^r 


zikl 


where the projection tensor A lk is defined as 


(50) 


A ik = g ik - U i U k . 


(51) 


The generalized 6-indices 8 —Kronecker tensor 8 l ^ ns (see, e.g., [T] ) may be defined by 
a recurrent formula through the 8 —Kronecker tensor with four indices, 8^ n , as 


8 ikl = 8 

mns ' 


8 kl +8 i 8 kl +8 l 8 

mns 1 nsm 1 


kl 

mn ? 


rik — ri rk _ ri r, 
u mn u ™ u -™ u ' 


5 k 

'm n 


k 

n^m ' 


(52) 
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Upon contraction, equation (EDI) yields another useful identity 


= S‘LU,U‘ = -A!„ . (53) 

The tensors £i k and (/i -1 )ifc are symmetric, but is in general non-symmetric. The 
dot denotes the position of the second index when lowered. These three tensors are 
orthogonal to U\ 


£ikU k = 0 , (fi~ 1 ) ik U k = 0 , v k U l = U = v k U k . ( 54 ) 

Using the equation (TIBI) , one can show through straightforward calculations that 

£ im = A im + 2 X ikmn U k U n , 

(ft- l )pq = A pq ~ lVpikX ikmn Vmnq = A pq ~ 2 *X* p l qs U l U S , 

v p m = VpikX ikmn U n = *X p in m U l U n , ( 55 ) 

which in turn satisfy the relations (USD- From the relations given in EHD one sees 
that the non-minimal interaction of the gravitational and electromagnetic fields 
effectively changes the dielectric and magnetic properties of the vacuum, and pro¬ 
duces a specific magnetoelectric interaction. In this sense, under the influence of 
non-minimal interactions the vacuum behaves as a material medium, called a quasi- 
medium. Note from EHD that the tensor x lkmn predetermines the changes in the 
dielectric properties of this quasi-medium, the double-dual tensor *x p i qs influences its 
magnetic properties, while the dual tensor *X P i n m produces magneto-electric effects. 

In order to complete this analogy, one can write the relationships between the 
four-vectors electric induction D l and magnetic field H\ on one hand, and the four- 
vectors electric field E l and the magnetic induction B l on the other hand. These 
relations are 


D i = e im E m - B l u, \ , Hi = v, m E m + (fi~ 1 ) im B m . (56) 

The vectors D l , H l , E l and B l are defined by the following formulae: 

D l = H ik U k , W = H* ik U k , E l = F ik U k , B i = F* ik U k . (57) 

These vectors are orthogonal to the velocity four-vector U l , 

D% = 0 = E% , H% = 0 = B% , (58) 

and form the basis for the F rnn and H mn tensors decomposition 

Fmn E m U n E n U m rjmniB , H mn D m U n D n U m Vmn iH . (59) 
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3.2.4 Master equations for the gravitational field 

We are working with a non-mininral electro-gravitational system, with the coupling 
terms linear in curvature, with the additional restrictions that the system is also 
linear in the Maxwell tensor, and the gravity part is Einstein-Hilbert. In this non- 
minimal theory, linear in the curvature terms, the equations for the gravity field (J 2 Hj) 
can be written in such a way as to look like the standard form of Einstein equation, 
i.e., as 

Rik - \r 9ik = « 4 eff) • ( 60 ) 

The effective stress-energy tensor T^ 1 in the right-hand-side of (1HUI) is quadratic in 
the Maxwell tensor and takes the following form 

bf 1 = bf+gib, 1 ’ +12 bf + is b? • (a) 

The linear part of the electromagnetic stress-energy tensor T® is given in equation 
da- The definitions for the other three parts of the stress-energy tensor, T^\ 
and T®, are 

= R T® - l -R lk F mn F mn - JfeV'VifF™/™) + k (F mn F mn ) , (62) 


rjl{2) 

^ ik 


2 9ik 


V V7 ( TTimn T7il \ 75 rr’inn rpl T? in f D P i r? 77 \ 

mV|(r -T.,J — ttlm-F r. n —r \iMirkn + J^kl-fin) ~ 


7imn rpl 


In ( 


-R mn F im F kn -^V l S7 l (F in F k n )+]^ l [Vi{F kn F ln )+\/ k (F in F ln )\ , (63) 


I O 1 

rn(3) _ 1 rymn.ls 77 77 ° jpls (Tp-n p i p-n p ^ 1 V7 V7 / r-n r-m 1 rp n T?-m\ 

1 ik ~ ^9ikM r mn ris--r (A Kknls + H in i s )-- V m V n (^i -\-r k t { )■ 

(64) 

Note that T ik in equation (1H3I) takes the same form as the stress-energy tensor 
constructed in [ 22 ] . 

While the stress-energy tensor of the electromagnetic field, T^\ has zero trace, 
the effective stress-energy tensor has a nonvanishing trace. Indeed, = 

g lk T^ s) , is given by 


J^(eff) _ _ 


9i 

-92 

-93 

1 


1 RF mnF ™ + lv k V k (F mn Fn 


R mn F k m F kn + -W k W k (F mn F mn ) + V m W n (F kn F krn ) 


R mnls F mn F ls + V m V„(F fcn F, 


kn : 


2 

mnls 


km ) 


2 X mnls FmnFis - (92 + 93)V m V n (F kn F km ) 


;(3gi + 92)V k V k(F mn F r ‘ 


(65) 
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Note that the sign of the trace is not defined a priori, depends on the specific model 
one uses. This feature also happens in non-linear electrodynamic models (see, e.g., 

E3). 

Equations m-m contain covariant derivatives of the Maxwell tensor only, and 
do not involve derivatives of the Riemann tensor, Ricci tensor and Ricci scalar. Thus 
for a given electromagnetic held they form a system of differential equations contain¬ 
ing second order partial derivatives in the metric. Nevertheless, these equations have 
to be completed by the self-consistent equations of non-minimal electrodynamics (©, 
(PI), and (HP, which contain the covariant derivatives of the Riemann tensor, Ricci 
tensor and Ricci scalar. In general, the Maxwell tensor, envisaged as a solution to 
equations (ED, and tup, depends on the second order partial derivatives of the 
metric. Thus, in general, the equations for the gravitational held become of fourth 
order. However, the parameters qi, q 2 and are arbitrary and may be hxed in an 
appropriate way. So, the question of whether or not there are models which are 
effectively of second order in the derivatives of the metric is pertinent. Below in 
section 3.3. we show explicitly one such a model. 

3.2.5 Bianchi identities 

Since the Einstein tensor in the left-hand-side of equation © is divergence-free, 
the effective stress-energy tensor m-m has to be conserved, i.e., 

v fc T M0 = 0 . ( 66 ) 

In order to check directly that this is true, one has to use, hrst, the Maxwell equations 
© and (1231) with (J2HD, and second, the Bianchi identities and the properties of the 
Riemann tensor, V iRklmn T iRikmn T ^kRlimn 0 and Rklmn T Rmkln + Rlmkn 0 , 
as well as the rules for the commutation of covariant derivatives, which for vectors 
yields (V,V fc - VfcV,)^ = W m R mlk . 

3.3 An example: static spherically symmetric gravitational 
and electromagnetic fields non-minimally coupled 

The line element for the static spherically symmetric model has the form 

ds 2 = B(r) c 2 dt 2 — A(r ) dr 2 — r 2 (d6 2 + sin 2 9 dtp 2 ) . (67) 

Assume also that the electromagnetic field inherits the static and spherical symme¬ 
tries. Then the electric held potential A has the form A = p{r)8^. The Maxwell 
tensor happens to be equal to F lk = p'[8\8\ — tf^,), where a prime denotes the 
derivative with respect to r. To characterize the electric held, it is convenient to 
introduce a new scalar quantity E(r) as follows, 

E 2 (r) = —E' l Ei = --F ik F ik = —F 2 = —p' 2 , ( 68 ) 

\ ) * 2 k AB r0 AB r ’ v ; 
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where the four-vector E l is defined in equation ED, and the velocity four-vector is 
chosen to be equal to U l = 5 l 0 B~z. 

To fix the sign we choose F r0 = —(AB)^E{r) and F r0 = (AB)~^ E(r). For 
this electromagnetic field the Maxwell equations © are satisfied identically, while 
equations ED and ED give only one non-trivial equation when i — 0, 


r 2 E(r) (l + 2x% r (r)) 


0 . 


(69) 


The function E{r ) can then be found to be 

E (r) = 2 ® ( t , where e r r (r) = 1 + 2y% ( r ), (70) 

r z £ r r {r) 

and Q is a constant. Assume now that the space-time with metric ED is asymptot¬ 
ically flat and y°Q r (oo) = 0. Then, the constant Q in ED coincides with the total 
charge of the object if (p(r) —> — at r —> oo. Using (ED one can compute the term 
X°or(n), 



(?i + Q2 + Qs) 


B" 

2 AB 


+(2?i + q 2 ) 


1 

2 rA 



(. B ') 2 A'B' 
AAB 2 ~ AA 2 B 



+ 



(71) 


So, from equation ED, one sees that generally, E(r) contains derivatives of the 
metric up to the second order. With such an electric field, equations ED - flHD for 
the gravitational field become of the fourth order. To illustrate this statement take 
the trace of equation ED, E = —^T^ eS \ where the trace T ( eff ) is given in ED- For 
the metric ED and the electric field ED-ED the trace equation takes the form 

1 \B" (. B') 2 A'B' 2 (B' A'\ 2 

k \ ~B 2 B 2 ~ 2AB + r\B~^)~^ 2 ^ [ ~ 


= (E 2 n3q 1+ 2q 2 +q 3 )+(E 2 )' 


(B' A 1 2 \ 2 

(3gi+2g 2 +g 3 ) 7T77“TT+ _ +-(^ 2 +^ 3 ) 


+E 2 


(gi+?2+g 3 ) - 


B" (B 
t v 


/\2 


2 B 2 A 
A'B' 2 N 

H-7T 


+ 


B 2 B 2 2 AB 




(72) 


Generally, equation ED includes the first and the second derivatives of the square of 
the electric field E(r ), which contains, in its turn, the first and the second derivatives 
of the metric coefficients. Thus, for generic qi, q 2 and q 3 we obtain a fourth order 
scalar equation for the gravity field. Direct calculations show that the equations 
derived from ED for the sets of indices tt, rr, 69, <p<p display the same features. 
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Now, when the susceptibility tensor is proportional to the double-dual Riemann 
tensor, i.e., q± + q 2 + (fe = 0 and 2 q 1 + q 2 = 0 or qi = g 3 = — q and q 2 = 2 q, all the 
derivatives disappear from the expression for E(r), providing the formula 


E{r) 


Q 

r 2 + 2q(l- i) 


(73) 


Thus, we recover the result obtained by Miiller-Hoissen and Sippel in |2£[ for the 
special model with qi = q 2 = 7, = —2y. Moreover, equations dnnp-dMD simplify 

significantly, in particular, equation (E2D yields 


R 


2 nq 
r 2 A 


r(E 2 )'+E 2 (2-A) + ^E 2 (R 



(74) 


This equation is, evidently, of second order with respect to the derivative d/dr. For 
/L(oo) = 1 this electric field is asymptotically Coulombian. Formally, (I7H1) has a 
form of the type discussed in mEmnnm. We intend to consider such a model in 
a future work. 


4 Conclusion 

We have established a new self-consistent system of equations for the gravitational 
and electromagnetic fields. The procedure was based on a non-minimal and non¬ 
linear extension of the standard Einstein-Hilbert—Maxwell Lagrangian. The class of 
systems we have studied includes non-minimal electro dynamic equations, containing 
the Riemann and Ricci tensors and the Ricci scalar both in the non-linear and linear 
versions. 

This class of models of non-minimal and non-linear coupling of the gravitational 
and electromagnetic fields is of great interest, since the appearance of cross-terms in 
the Lagrangian leads to modifications of the coefficients involving the higher-order 
derivatives both in the Maxwell and Einstein equations. This means, in particular, 
that the velocity of the coupled gravito-electromagnetic waves should differ from the 
speed of light in vacuum. 

The general field equations obtained in the paper can in principle be classified 
using the explicit dependence of £(/( a &)) on J( afe ) in the non-linear theory, whereas 
in the linear theory one uses the phenomenological parameters q \, q 2 and q%. This 
is important for two reasons. First, one should search for non-minimal models 
in which the gravitational field is described by equations of the second order in 
the derivatives of the metric. We have shown explicitly, that static spherically 
symmetric configurations satisfy such a requirement if the susceptibility tensor is 
proportional to the double-dual Riemann tensor. This model requires a detailed 
analysis and we intend to consider it in a separate paper. Second, for the non- 
minimal non-linear coupling between electrodynamics and gravitation one should 
search for master equations (no matter whether they are of second or of higher order) 
admitting non-singular, regular, solutions for the gravitational and electromagnetic 
fields. 
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